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Abstract In a previous paper (S. Ghosal and Z. Chen Bull. Math. Biol. 2010 72, 

pg. 2047) it was shown that the evolution of the solute concentration in capillary 

electrophoresis is described by a nonlinear wave equation that reduced to Burger's 

equation if the nonlinearity was weak. It was assumed that only strong electrolytes 

(fully dissociated) were present. In the present paper it is shown that the same 

rG • governing equation also describes the situation where the electrolytic buffer con- 

Mh' sists of a single weak acid (or base). A simple approximate formula is derived 

Q . for the dimensionless peak variance which is shown to agree well with published 

• ^ ' experimental data. 

C/3 , Keywords capillary electrophoresis • electromigration dispersion 

o 

*'0' 1 Introduction 

Capillary electrophoresis (CE) is a technique for separating a mixture of macro- 
ions in aqueous solution by exploiting the fact that the migration velocity in an 
, applied field depend on the size and charge of the molecule. It is a widely used 

J^ I laboratory tool in bio-analytical chemistry. Further background information may 

vQ , be found in the authors' earlier paper in this journal [1] (henceforth refereed to 

f^ ' as GC) and some of the references cited there. There are also several textbooks 

devoted to the subject ^2ll3j. 
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The mathematical formulation of the problem consists of a set of coupled equa- 
tions describing the transport of ions in response to electric fields and diffusive 
fluxes. The electric field in turn is determined by the concentration distribution 
of ions. Coupling to hydrodynamics could occur due to electro-osmosis, but is ne- 
glected here for simplicity. Electromigration dispersion (EMD) is caused by vari- 
ations in the local electrical conductivity in the vicinity of the solute peak which 
gives rise to a nonlinearity in the equation for solute concentration. 

Simple one dimensional mathematical models of electromigration dispersion 
neglecting the effects of diffusion have been studied by various authors [ll[5l|6]. 
The problem is reduced to a single nonlinear hyperbolic equation for the concen- 
tration of sample ions. Analytical and numerical solutions describe the character- 
istic wedge shaped profile observed in experiments [71|8]. The restriction to zero 
diffusivity was removed recently by the authors [T] who considered a 'minimal 
system' of three ions - the sample ion, co-ion and counter-ion - all being strong 
electrolytes (fully dissociated) . The diffusivities of the three ionic species were con- 
sidered non zero but equal. The sample concentration was then shown to obey a 
one dimensional nonlinear advection diffusion equation which reduced to Burgers' 
equation if the sample concentration was not too high. Thus, the concentration 
profile could be obtained analytically as a function of position and time, providing 
useful insights into the nature of electromigration dispersion. 

The solutes of interest in CE are often biological molecules for which the charge 
is quite sensitive to the pH of the surrounding electrolyte. Many molecules are also 
unstable outside a narrow pH range. Thus, CE must be performed in a medium 
where the pH is kept as constant as possible. In order to maintain a stable pH 
an electrolytic buffer containing a weak acid or base is used as the background 
electrolyte. The buffer often contains several ionic species as well as other additives 
to achieve different functions (e.g. prevent adsorption to capillary walls). Thus, 
the simplifying assumption made by GC that the buffer consists of a single strong 
electrolyte is often not consistent with laboratory practice. 

In this paper we show that the one dimensional model derived by GC for strong 
electrolytes may also be applied to an idealized model of a buffer consisting of a 
single weak electrolyte. The only change required is in the definition of a parameter 
in the model that characterizes the strength of EMD effects. A simple formula for 
the "number of theoretical plates" - a dimensionless measure of peak dispersion - 
is derived and is shown to agree well with published experimental data. 

The rest of the paper is organized as follows; the theory of EMD in the presence 
of a weak electrolytic buffer is provided in the next section followed by a discussion 
of the experimental work and comparison with the theory. A summary of our 
results and discussion of the validity of our underlying assumptions is provided in 
the concluding section. 



2 Theory 

For definiteness, we consider an acidic analyte, HnA + nH20 — >■ nH^O^ + A"" in 
aqueous solution buffered with a weak acid HX + H2O ^ H^O'^ + X~ . The auto 
ionization of water is neglected. The analyte is considered to have a fixed charge 
{ze), but, the buffer could exist in either a charged {X~) or neutral {HX) state. 
Cationic analytes or buffers with multivalent ions can be accommodated easily in 
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our analysis but we only discuss the simplest situation for convenience. Further, we 
assume that all ions have the same mobility (u), and therefore, identical diffusivity 
(D), and that the transport problem is entirely one dimensional. In particular, we 
assume that electroosmotic flow, when present, may be described as an advection 
of the ions with a constant velocity mq in the axial direction. We will return to the 
question of the validity of these assumptions in the concluding section. 



2.1 Derivation of a reduced system 

The coupled equations describing the concentrations of hydrogen ions (c-|-), acid 
ions (c_), sample ions (c) and the neutral form HX (cq), are then 

^ + |.[(.o + e.i.)c,]=i.^+., (1) 

^ + |.[(.o-e.i.)c_]=i5^+., (2) 

P,+l[i-o + .eum=D^, (4) 

where E is the local electric field, x is the distance along the capillary, t is time, e 
is the electronic charge and r = k^co — kac+c- is the net dissociation rate of HX. 
Here, ka, k^ are constants that characterize the rates for the forward and reverse 
reactions of the weak acid. We will assume that the x-axis points in the direction 
of peak motion from the inlet towards the detector. These equations differ from 
those considered by GC in the presence of the source terms on the right hand 
sides of the first three transport equations and in the existence of the additional 
variable, cq. Since characteristic spatial scales are always much larger than the 
Debye length, local electro-neutrality holds. Thus, 

cj^ — c- -\- zc = 0. (5) 

If we multiply equations (Il])-(I3|) by the respective ionic charges, sum them, and use 
equation ([5]) , we get an equation that describes the constancy of electric current 

^ [e^u{c+ + c_ + z^c)E\ = 0. (6) 

Note that the net contribution from the diffusive fluxes as well as those from the 
source terms vanish exactly. Equation (|6]) may then be integrated; 

{c+ + c-+z'^c)E = 2cocEoo, (7) 

where Eoo and Coo are the electric field and cation (or anion) concentration re- 
spectively far away from the peak. 

Equations Q and ([7]) provide two algebraic relations among the five depen- 
dent variables cj^,c^,co,c and E. To reduce the system of equations (HI - (HI) and 
([6]) to a single one dimensional equation we must seek more such algebraic rela- 
tions. A third relation is provided by an approximation first introduced by Saville 
and Palusinski [9': the time scale associated with the dissociation-recombination 
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reactions represented by the last term on the right hand side of equation (JS)) is so 
small compared to all transport time scales that equation ^ niay effectively be 
replaced by 

r = k^co — kaC+c- = (8) 

or kco = c-\-C-, where k = kj^/ka is the acid dissociation constant. Thus, locally, 
the acid is in equilibrium with its dissociation products. 

In GC it was shown that the Kohlrausch function K = {c+ + c_ + c)/u is a 
passive scalar that spreads only by diffusion (or, if electroosmosis is present, is 
also advected with a constant velocity uq). Therefore, electrophoretic migration 
relative to the fluid rapidly advects the solute peak into a region where K is 
essentially equal to its far field unperturbed value. This provided the additional 
algebraic relation (the constancy of K) that finally enabled the reduction of the 
transport equations to a one dimensional system. This approach does not work in 
the current problem, as it may be readily shown from the transport equations ([T|). 
^ and (|3|), that the evolution equation for K, 

^+uo— = D— + — (9) 

dt dx dx^ u 

now has a source term, and K is therefore no longer a passive scalar. 

Fortunately however, a fourth algebraic relation is obtained if one assumes that 
the buffering action of the weak acid HX is "perfect" so that 

C+ = Coo. (10) 

The justification of the approximation (|10p is rooted in the theory [TO] of acid- 
base equilibria where it is shown that in a mixture of a strong {HA) and weak 
{HX) acid, the perturbation in the hydrogen ion concentration is small when the 
buffering capacity of the background electrolyte is large. The buffering capacity 
is a maximum, if the solution pH (-logc+) is equal to the pKa (-logfc) of the 
acid and drops sharply if the pH differs by more than one unit from the pKa . In 
practice, the buffer consists of a mixture of the weak acid and its salt so that the 
pH and buffer concentration can be independently controlled. At a given pH, the 
buffering capacity increases with buffer concentration. The buffer concentration 
cannot be made too large however, as this increases the solution conductivity 
thereby reducing separation efficiency. In the context of our model represented by 
equations ([1]) - (|1|) and ([6]) , the significance of the assumption (|10p becomes clear 
if we use equations ([5]) and ([8]) to express c-|- in terms of the concentrations c and 
Co: 

c+ = v^feco + z2c2/4 - zc/2. (11) 

Equation (|lip may be replaced with equation (llOf) . if 

kco » z^<?, (12) 

that is, if the amount of undissociated acid in the sample zone is high enough to 
act as a "buffer" against pH variations. 

With four algebraic relations among the five dependent variables, equations ([1]) 
— ((31) and (O may now be reduced to a single equation for the normalized solute 
concentration (f> = c/coo : 



dt dx 



UQ + 



1 — a(f> 



^S. (13) 
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Fig. 1 Comparison of the normalized effective diffusivity given by equation 1151 1 [solid line] 
writh that predicted by the approximate form equation II17II [broken line]. 



where v = zeuEoo is the electromigration velocity of an isolated solute ion and 

a = -i2(l + 2). (14) 

Equation (|13(1 is identical to the corresponding equation for strong electrolytes 
(equation (13) in GC; except there mq =0). The only difference is in the definition 
of a which is given by equation (|14p instead of equation (14) of GC, which is 
a = (1 — z'^)/2 for a univalent background electrolyte. 



2.2 Dispersion 



Equation psp is analytically solvable in certain limits and its consequences were 
discussed at length in GC. Here we are interested in the peak dispersion of the 
solute as measured by A*' = L^ja^ where L,^ is the distance between the injection 
point and the detector, and, a^ is the variance of the concentration when the peak 
reaches the detector. One of the consequences of equation (|13p is that if <^ is not 
too large, the variance of the peak increases in proportion to the time so that the 
spreading may be described in terms of an effective diffusivity 

°eff-2I.(|-|) (15) 

where Fn is the n th moment of a certain function F(x) defined by equation (23) in 
GC. The degree of sample loading is characterized by a length r = J_°° ^(a;,*) dx 
and the relative importance of diffusion is characterized by a Peclet number P = 
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r\v\/ D. The effective diffusivity fias the following asymptotic forms at small and 
large values of the Peclet number P : 

(D ifP«l, 

^eff~|i|«^|rifP»l. (16) 



For the purpose of practical applications, equation psp . which requires the nu- 
merical evaluation of integrals may be replaced by the simpler formula 

^efr = -^ + ^l"H^- (17) 

The ratio D^g/D given by equation (|15p depends solely on |q:|P. The approxi- 
mation to Da given by equation (|17p approaches the theoretical value given by 
equation psp in the asymptotic limits of small as well as large Peclet numbers. 
In Figure 1 we plot D^g/D as a function of \a\P evaluated using equation (|15p 
as well as the approximate form equation p7p . It is seen that the approximate 
form provides an excellent approximation for all values of |ci|P while having the 
advantage of algebraic simplicity. 

The number of theoretical plates as a function of the magnitude of the applied 
voltage, V, may be calculated from the classical result [2J 

N=iE^±^, (18) 

where ^o is the electroosmotic mobility and /i = zeu is the electrophoretic mobility 
of the sample ion. Equation psp may be generalized by replacing the diffusivity, 
D, by the effective diffusivity -Dgff- We also allow for an initial peak variance {a^), 
and, an inlet to detector distance {fL) that is slightly smaller than the capillary 
length (L). Thus, 

where the constants A and B are respectively 

A=^A^ (20) 

and 

We will now compare equation (|19p for A'' with published experimental data. 

3 Experiments 

Jorgenson and Lukacs [12_ derived equation psp as a simple model of peak dis- 
persion. The experimental data that they presented however, showed the expected 
linear dependence of A^ on V, but only over a limited range of V. At large V, the 
linear dependence was found to exhibit a saturation effect so that A' approached a 
limiting value rather than increase indefinitely. Jorgenson and Lukacs attributed 
this to the effect of Joule heating in the capillary. Delinger et al. |T3] later re- 
interpreted Jorgenson and Lukacs data and suggested that the saturation was 
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Fig. 2 Symbols are data from Figure 3 of Lukacs and Jorgenson |11| re-plotted as A'^ ^ vs. 
c* on the right panel. The solid line is the best fit linear regression N~^ = o, + 6c». 



more likely due to the initial variance {(jq) of the injected zone. They also pre- 
sented additional experimental data to support their interpretation. Equation (|2ip 
shows that electromigration dispersion is another possible reason for the observed 
saturation of A'' at large V . Whether the initial variance or the electromigration 
effects dominate will depend on the experimental parameters. As the voltage is in- 
creased further, Joule heating results in a decrease of A'^. In some experiments [Tl] 
Joule heating masks any effect of electromigration and no distinct plateau region 
is observable in the N vs. V curve. Instead, the A'' vs. V curve displays a single 
maximum. Thus, in any comparison of equation (|19p with experimental data, care 
must be taken to distinguish electromigration effects from those due to the initial 
variance and Joule heat. Equation H2ip suggests that this can be done in the clear- 
est manner if one varies the sample loading F while keeping all other experimental 
parameters invariant. In this situation a linear dependence of A'^^ with F should 
be observed. Such an experiment was published by Lukacs and Jorgenson [,11._ and 
we shall compare equation (|19p with their data. 

Specifically, we use the data shown in Figure 3 of the paper by Lukacs and Jor- 
genson [If (henceforth referred to as LJ) where the solute was Dansyl-isoleucine 
in a 0.05 M phosphate buffer at pH 6.86. In the experimental setup, the solute 
concentration (c*) alone was varied. In each run, the same length of plug {tj 
was injected into the capillary by electrokinetic injection. Thus, the initial profile 
4>{x,Q) was a square wave, and therefore, F = CtHjCoo and ctq = (^ jVl. Accord- 
ing to equation (|19() . A''"^ should have a linear dependence on the concentration: 
7V^^ = a-\- fee* where 



2D 



+ ■ 



1 



2, , Iwl (. 1 



9 fio + /J. f L Coo 



(22) 
(23) 
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Table 1 


Experimental parameters 


corresponding to the data 


shown in Figure [2] 




L (cm) 


f 


V (kV) 


/^o (cm^/Vs) 


/i (cm2/Vs) 


D (cm2/s) 


z 


Coo(M) 


100 


0.80 


30 


4.79 


X 10-4 


-1.37 X 10-* 


5.35 X lO-f' 


-0.65 


0.025 



The symbols in Figure [2] are the experimental data which does show that N ^ has 
a linear dependence on c*. Fitting a straight line through the data points gives 
a = 2.57 X 10"*^, b = 2.98 x 10"^ M"\ 

The reported values of the various constants appearing in equations (|22ll and 
(|23p are summarized in Table 1. The value of the electroosmotic mobility, ^q at 
pH 6.86 for pyrex was read from Figure 4 of LJ. The electrophoretic mobility, /i, 
diffusivity, D, and valencqj z of Dansyl-isoleucine reported in Table 1 were taken 
from the paper by Walbroehl and Jorgenson [15_, where these parameters were 
measured under settings identical to those employed in LJ. The buffer strength of 
our idealized buffer HX when pH = pKa, is, by the Henderson-Hasselbalch relation 
[HX] + [X^] = 2coo. Equating this to the reported buffer strength of 0.05 M gives 
the value of Ccx3 indicated in Table 1. Unfortunately, the injection Voltage (V*) and 
injection time (r*) for the electrokinetic injection process were not reported, so 
that, £ is not known. However, if we use the value of a obtained from the best fit 
line in Figure[2] we get £ = 0.31 cm. From the relation I = {fj, + ho)V*t^,/L we can 
determine that V*r» ^ 90 kV • s, which appears consistent with what is typically 
reported in similar experiments [13] . If we substitute this value of the injection plug 
length, £, into equation (|23|) and use the value of b found from the linear regression, 
we deduce that a = 0.22. The phosphate buffer is more complex than the idealized 
model of a buffer considered here. However, if we put z = —0.65 in equation p4p 
we get a Ri 0.1. If the background species was a fully dissociated 1-1 electrolyte we 
would have (GC) a = (1 — a^)/2 f» 0.3. These numbers are commensurate with the 
value computed from the experimental data using equation (|23p . which suggests 
that the value of a is probably not very sensitive to the simplifications adopted to 
model the background electrolyte. 



4 Conclusion 

The problem of electromigration dispersion of a solute was considered in the pres- 
ence of an idealized 1-1 weak electrolytic buffer. The analysis complements earlier 
work (GC) where the background electrolyte was regarded as fully dissociated. It 
was shown that in both models the solute transport is described by the same one 
dimensional transport equation; they differ only in the definition of the parameter 
a characterizing the strength of EMD effects. A simple expression for the number 
of theoretical plates was derived by replacing the analytical expression for the ef- 
fective solute diffusivity by an approximate form. This expression for A'^ was then 
compared with published experimental data and good agreement was found. 

The observed agreement between theory and experiment requires further ex- 
planation, since the phosphate buffer used in the experiment is a complex multi-ion 



^ Dansyl-isoleucine has a single negative charge on the carboxyl group but its effective charge 
is likely reduced by shielding effects. 
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mixture and does not correspond to either the strong electrolyte model of GC or 
the idealized weak electrolyte model considered here. The explanation must be, 
that, the fact that the concentration is governed by equation (|13p is independent 
of the specific model assumed for the buffer. This is quite plausible, since to derive 
equation (|13|) all that is required is that one must be able to find A'' — 1 linear 
algebraic relations among the A'^ chemical species in solution. The analysis pre- 
sented in this paper and in GC are simply two alternate ways in which this can be 
accomplished through different assumptions about the nature of the background 
electrolyte. There are of course other ways to a similar end. For example, in a 
more complicated multi-component background electrolyte, in addition to local 
electro-neutrality, one may assume local equilibrium for a subset of the species 
and set the concentrations of certain other species to zero or to a constant value 
depending on the details of the buffer chemistry. Alternatively, one could take a 
more heuristic view point and simply say that the local migration speed of a solute 
ion is a function of the local solute concentration, and, if the solute concentration 
is not too high, this function may be linearized: Ve{cj)) = ^^il + '^4' + ■■■]• In this 
approach a would be an empirical parameter. Our analysis suggests, that in gen- 
eral, a may be well approximated by a parabolic function of the solute valence 
a = ao -\- ciiz -\- a^z where the coefficients aQ,a\^a2 depend on the buffer com- 
position. For a 1-1 strong electrolyte, we found (GC), zq = —22 = 1/2 and z\ = 0, 
whereas for a 1-1 weak electrolyte, zi = 2:2 = —1/2 and zq = 0. Determining the 
coefHcients ao, Q^i a-nd a-2 experimentally may be an efficient way of characterizing 
the dispersive properties of electrolytic buffers in the laboratory. 



In arriving at equation (|19(1 . the local ion migration velocity was linearized: 
We(<^) = ^'/(l — oi(^~) ~ v[l + a(j)\ by neglecting quadratic and higher powers of 
(f). This is a valid approximation if <C 1. Furthermore, the smallness of is 
also inherent in the assumption of a "perfect buffer", equation pop . Indeed, if we 
require that equations pop . ([5]) and ([8} be consistent with the requirement p2p . 
then we must have 1 + 20 3> z 4> , that is, (l> must be small in comparison to unity. 
In the experiment, even at the highest solute concentrations used, (f) ~ c»/coo = 
0.006A//0.025M = 0.24. Thus, lyi <C 1 is a valid approximation in the experiments 
of LJ. 

In the current analysis as well as in GC, the diffusivity of all ionic species 
were presumed equal. This is not an entirely unreasonable assumption, since the 
diffusivities decrease weakly with molecular mass. For example, the diffusivity of 
Dansyl-isoleucine is only about three times larger than that of the sodium ion 
though its molecular mass is about sixteen times as great. The exception is the 
hydrogen ion which has a diffusivity about a factor of ten higher than Dansyl- 
isoleucine. However, in the phosphate buffer the cations responsible for conducting 
current are primarily the sodium ions, the hydrogen ion concentration is actually 
very small. For macro-ions with molecular weights in the range of kilo Daltons as 
well as for colloidal particles, differential diffusivity could turn out to be important. 
Such effects can be accounted for in our theory by including the diffusive current in 
the equation for current conservation, equation Q. This would make the electric 
field dependent not only on the solute concentration but also on its gradient. The 
resulting one dimensional model would include the effects of differential diffusion, 
but one may lose the convenience of having an analytically solvable governing 
equation. 
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The pyrex capillary employed in the experiment generated strong electroos- 
motic flow which was accounted for in a simplistic manner as a uniform advection 
along the capillary. However, the perturbation in the electric field generated by the 
solute peak also perturbs the slip velocity at the capillary wall which results in ra- 
dial shear and consequent Taylor- Aris dispersion by a well known mechanism [1^ . 
It may be shown [IT] that this effect can be accounted for by adding a contribution 
to the diffusivity that is quadratic in 4>. In the present experiments, this effect is 
therefore expected to be negligible due to the smallness of 4>- 

We have shown that the theoretical framework developed here and in our ear- 
lier paper (GC) is useful for analyzing real experimental data, not withstanding 
the fact that the theory pertains to very idealized situations. In particular, equa- 
tion (|19[) is a simple generalization of the formula (|18p introduced by Jorgenson 
and Lukacs p2| that has proved to be extremely useful in laboratory practice. It 
is hoped that equation (|19(1 could be used in a similar manner to characterize the 
effects of EMD. 
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